Traveling wave solution of the Reggeon Field Theory 

Robi PeschanskEI 

Institut de Physique Theorique, 
CEA, IPhT, F-91191 Gif-sur-Yvette, France 
CNRS, URA 2306 



We identify the nonlinear evolution equation in impact-parameter space for the "Supercritical 
Pomeron" in Reggeon Field Theory as a 2-dimensional stochastic Fisher-Kolmogorov-Petrovski- 
Piscounov equation. It exactly preserves unitarity and leads in its radial form to an high energy 
traveling wave solution corresponding to an "universal" behaviour of the impact-parameter front 
profile of the elastic amplitude; Its rapidity dependence and form depend only on one parameter, 
. . . the noise strength, independently of the initial conditions and of the non-linear terms restoring 

0^ ' unitarity. Theoretical predictions are presented for the three typical distinct regimes corresponding 

' to zero, weak and strong noise. 

o : 

CN . PACS numbers: 

■ I. INTRODUCTION 

(N : 

^ \ The question of the high-energy behavior of soft hadron-hadron amplitudes and in particular of their expanding 
i impact-parameter disk with rapidity is a rather old subject, but still not solved, being in the basically unknown realm 
Q-i' of non-perturbative QCD. However, one promising theoretical approach at early time is the Reggeon Calculus [H and, 
yi^ , in the formalism which we will be dealing with in the present work the Reggeon Field Theory (RFT) 0, where the 
^ ' amplitude is described in terms of an effective quantum field theory of "Pomeron fields" . It derives from a Lagrangean 

■ in 2H-1 dimensions, where "space" is the 2-d impact-parameter h and "time", the overall rapidity Y. In the studie 
performed during the 70's, and after a series of works dedicated to the renormalization group approach to the RFT 
Q , it appeared also d, 0, |^ that a physically interesting case is when one considers instead a "supercritical" bare 

^ ■ Pomeron P, i.e. when the intercept a is greater than 1 (in fact a > Oc > 1, where Uc includes the quantum effects of 
Cn the renormalizable RFT when the Pomeron field is at criticality 0)- In that case, unitarity is violated by the bare 
Pomeron (equivalent to a Born term) but expected to be recovered thanks to damping Pomeron interactions. It was 
shown that the impact-parameter disk was expanding like the rapidity Y, expressing a dynamical instability ofthe 
: RFT 5]. However, the field theoretical techniques known at that time did not seem to give much more indication on 
the solutions. 

, The notion of a Reggeon Field Theory and its use to describe "soft" hadronic interactions appeared, since that 
time and until recently, under various forms which may differ from the original version depicted in 0. For instance, 
they may refer to the similar approach deduced or inspired by QCD and its dipole model ^, 't']. Also, there exists a 
phenomenological interest for using a "supercritical Pomeron" in models based on interacting Pomerons Q. So, we 
want to specify now in which sense we use the RFT and what is different in our approach from the previous ones. 
First, we want to address the problem of finding the solutions of the full 2-dimensional transverse space problem 
including an explicit form of the impact-parameter dependence of the elastic amplitude. To our knowledge, explicit 
solutions have been only found in the zero-dimensional approximation only. Two-dimensional formulations of the 
QCD Pomeron calculus in the dipole approach are also widely discussed 0,0], but explicit solutions seem to be 
difficult to acquire. So, we restrict our analysis to the initial formulation |2| and thus our starting point for the 
Lagrangean is the original one 2]. Hence, when we use the term Reggeon Field Theory (RFT) we refer in the present 
paper to that precise formulation, up to a suitable generalization to be discussed further on. 

The goal of our paper is to update the original study of RFT by introducing new powerful tools known under the 
name of "traveling wave solutions" of non-linear evolution equations and already used in a different context for QCD 
evolution equations [9] . In fact one of our motivations is to try and give a theoretical answer to an old question raised 
by the phenomenology of elastic hadronic reactions, and its approach by interacting supercritical Pomerons. The 
"soft" hadronic elastic amplitudes seem to follow a common behaviour at high-energy, independent of the reaction 
one considers. While this property could be understood by the factorization properties of a single Regge pole exchange 
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(see, e.g. [10|), it is not known how a property may be obtained from an interacting Pomeron framework where the 
"bare" Pomeron input is deeply modified by the interactions. 

The relation of RFT with non-linear evolution equations (which will allow to use the traveling wave framework) 
has appeared since long in relation with statistical mechanics of out-of equilibrium processes. Starting with the deep 
relation between the original RFT with directed percolation there is a long list of works using this kind of 
connection , in particular with stochastic evolution equation of Langevin type. In fact some of the works (see, e.g. 
0,01) are using this connection to try and derive the stochastic evolution equation corresponding to QCD Reggeon 
Calculus. We shall indeed use some tools from statistical mechanics, in particular those developed in Ref. 12] to 
transform the RFT formulation in terms of a Langevin equation of known type and possessing traveling wave solutions. 

The main feature of traveling wave solutions of non-linear evolution equations is that they lead to "universality" 
properties, that are properties which will be of general value, i.e. irrespective of particular initial conditions or on 
features of the equation such as the form of the non-linear damping terms responsible for the unitarity restoration. 
Our hope is thus to provide through the traveling wave approach, an explicit high-energy solution of the RFT and in 
the same footing a physical understanding of the empirical "universal" properties of soft scattering amplitudes at high 
energies which are difficult to explain in a supercritical Pomeron framework. Our paper thus contains the theoretical 
derivation of the solutions of the (2-f 1) dimensional RFT [2] via the identification of the related stochastic non-linear 
Langevin equations. 

The paper is organized as follows. In section|TTl we show that the RFT, with a supercritical bare Pomeron as input, 
can be found equivalently realized by the 2-dimensional version of the stochastic Fisher and Kolmogorov, Piscounov, 
Petrovsky (sFKPP) equation, and that the elastic amplitude is solution of its reduction to the 1-dimensional radial 
(azimuthally symmetric) form. In section [Till we derive the main feature of the mean field (or deterministic) radial 
FKPP equation: the existence and universal properties of circular traveling wave asymptotic solutions. In section 
IIVI we introduce the effect of stochasticity by analyzing the solution dependence on the noise term. It appears with 
two markedly different regimes at weak and strong noise strengths. In section IVTl we present our conclusions and an 
outlook on the theoretical implications of the traveling wave picture. In the appendices we show the derivation of the 
solution in the deterministic radial case and an overview on possible phenomenological implications. 



II. FROM REGGEON FIELD THEORY TO THE 2-D SFKPP EQUATION 

The RFT with a supercritical Pomeron is defined 2] from the following ingredients, namely one propagator V ^ V, 
with coupling fi corresponding to the bare supercritical Pomeron intercept 1 + fj, > 1. There is a kinetic term in 
impact-parameter space with coupling identified with a', the slope of the bare Pomeron trajectory. For the Pomeron 
interaction vertices, one includes the triple Pomeron vertex, which gives rise to two possible contributions, i.e. the 
merging triple Reggeon term V+V — *■ V and the splitting term V — » V+V, with initially equal strength A corresponding 
to the triple-Pomeron coupling. 

The field theory action is defined in terms of quantum bosonic fields f and the conjugate (p by the action Q 

S[ip,^]^^ J(fb dY {^[dY-aV^]ip-fiipip-i\{(f(p^+^^ip)} . (1) 

As discussed in the imaginary coupling constant makes this theory non-hermitian and thus the fields (p and (f 
do not play a symmetric role through time reversal. It is indeed convenient to perform the field transformation 
iip, (p — !■ —i(p, giving rise to the modified action 

S[ip^(p] = ^ Jd^b dY {>p[dY~a'S/^]ip-fi ^ip + Xipip^- X^'^ip} . (2) 

One important outcome of RFT is its remarkable connection with problems of non-equilibrium statistical physics [ll[ . 
Following a known technique [T3|, the fields (p can be integrated out, and play the role of auxiliary fields appearing 
as external source fields for the deterministic part (for linear terms in (^)and the noise terms (for quadratic terms 
in If) of a non-linear Langevin equation, as we shall now see. Following |12l |. one linearizes the remaining quadratic 
(p'^ contribution in Q by introducing a stochastic white noise via a Stratonovitch transformation^, in such a way 
that all terms become linear in ip. Then performing the path integral over ip boils down to a a nonlinear Langevin 



^ The linearization of the quadratic terms in the action is a well-known procedure. The interested reader will find in rcf. '12'| a detailed 
derivation of the transformation of the action l(2j leading to the Langevin formulation (3i . 



3 



equation for the field Lp which now acquires the interpretation of random reahzations of the (properly normalized) 
elastic scattering amplitude T, namely 

^ T (y, hj = a' VlT + nT-\T^ + V2a'X T v{Y, 6 ) , (3) 
where the white noise verifies 

iv{Y., b ), iy{Y', b' )^ = S{Y' - Y)5'^{b' ~h) . (4) 

We have now to introduce the appropriate normalization of the scattering amplitude in impact-parameter space 
which is imposed by the unitarity limit T = 1. This comes as a constraint both on the deterministic and on the 
stochastic part of ([3]), since T = 1 should appear as a stable fixed point of the equation^ (the other fixed point is 
the "unstable" fixed point at T = 0, since the rapidity evolution increases, at least in average, the value of T). The 
unitarity constraint thus leads us to modify Eq.Q to get the following form 

^T{Y,b) =a' VlT + M (T - T^) + ^2a'Kii (T - T^) v{Y, b ) , (5) 

where it imposes equal coupling A = to the terms in T and of the deterministic equation and adding a term 
in the noise factor in order to ensure it to vanish at the black disk limit. 

A key feature of Eq. ^ compared to the initial formulation (j3|) is the introduction of the parameter k, which plays a 
crucial role both physically and mathematically on the nature of the solutions. At first we note that unitarity imposes 
no a priori constraint on the noise strength and thus allows for the introduction of the parameter k. Physically, k 
introduces a parametric factor between the strength of the merging term V+V V and the splitting term P V+V. 
This degree of freedom may come from at least two physial motivations. First, we will see that the traveling wave 
solutions possess "universal" features, in particular they will remain valid for more complicated merging factors 
{e.g. T"',n > 2 or even with a positive monotonous function T'^f{T) with /(I) = 1). Hence there is a priori no 
constraint of equal coupling between merging and splitting terms. A more intringuing motivation may come from 
an analogy with the dipole picture. Taking into account their size, the merging and splitting rules for dipoles imply 
a size-dependence of their effective coupling strength. Indeed, "fat" dipoles may merge more easily than "thin" ones, 
since it requires a matching of their transverse coordinates, while splitting does not seem to require such an effect. All 
in all, we find it physically suitable to consider the generalized Eq.® as the basic equation to be solved. Obviously 
taking k = 1, we recover the original RFT, up to a quadratic term in the noise which is easy to reinterpret as a 
four-vertex, see further. 

Our basic starting point is to point out that Eq.([5]) is (by introducing canonical variables, see further) the extension 
in two spatial dimensions of the Fisher and Kolmogorov, Piscounov, Petrovsky (FKPP) equations (for k = 0), or (for 
K ^ 0) its stochastic extension (sFKPP), see 

Refs.|i,[il,[il[il[i3,l23- 

This will allow us to find the solutions of the 
RFT for a supercritical bare Pomeron, thanks to modern tools'^ applied to the old and yet unsolved RFT problem. 

Mathematically speaking, we are looking for solutions which are not dependent of the forms of the non-linear terms 
provide they ensure a stable fixed point, that is T = 1 in our case. Universality also means that the solution is 
independent of the initial conditions after some "time" (here, rapidity) evolution interval. This defines a "universality 
class" of solutions, which will ultimately depend only of the value of k, that is on the noise strength. If we were in 
the situation of statistical physics at equilibrium , we could consider k as the order parameter of the problem. In our 
case it will allow to separate different regimes (but not necessarily separated by critical points.) We also note that the 
quadratic term in the noise can be easily reinterpreted in the field theoretical framework as a V+V — *■ V+V coupling 
in the RFT framework. This is equivalent to the following RFT action 

S[ip,^] = -^ Jd^b dY {^[dY-a'\7^]ip-fi i^{ip- ip^)-Kfi >f^{ip- (p^)} . (6) 

To complete the theoretical preliminaries it is worth mentioning that from the point of view of statistical physics, it 
is known that more general Langevin equations can in turn be analyzed in terms of a bosonic quantum field theory^ 



^ The stable fixed point is finally reached at infinite rapidity. In fact, one could technically equivalently consider an unitarity-preserving 
fixed point of Q at T = fi/X < 1. However, it is more often considered that the black disk limit T = 1 is the physical one. 

^ It is to be mentionned that a first connection between Reggeon Field Theory and circular traveling waves applied to cluster growth 
appeared already in 

Doi and Peliti, Refs. \2^, addressed the related problem of mapping master equations for reaction-diffusion processes to a field theory 
action. 
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FIG. 1: Traveling waves in 1 dimension. The traveling waves are asymptotic solutions u{x,i) of the 1-dimensional FKPP 
equation, translation invariant in time and joining the unstable fixed point [transparency i.e. dilute medium), to the stable 
fixed point [strong absorption i.e. dense medium). The figure is from Ref. [9|. 



[2l|. This formalism is particularly convenient to treat fluctuations superimposed to mean- field equations. Hence 
both techniques coming from statistical and particle physics can be joined together to get a deeper understanding of 
the original RFT problem and find the structure of its solutions, i.e. identifying its "universality class". 

Let us introduce now canonical variables allowing to put ^ in the generic form of the sFKPP equation. By suitable 
redefinitions 



T{Y,h)= U{t,f) ; ^i{Y-Yo)=t 




Eq. ([5]) can be recast in the canonical form 

^ U{t,r)^VlU+ U- U^ + e ^Uil-U) v{t,T) , (8) 

where the only remaining dimensionless parameter defines the normalized noise strength e as a function of the product 
of the "splitting over merging" factor k, and the "super-criticality parameter" /i. Eq. ([5]) is the canonical form of 
the nonlinear stochastic Fisher-Kolmogorov-Petrovski-Piscounov (FKPP) equation. It is worthwhile to note that the 
time relation in ((7|) is defined up to a rapidity translation Y —Yq. 

To remind known properties in dimension one, the remarkable feature of the FKPP class of equations [ol. [T3. flTj 
is to admit asymptotic traveling wave solutions, i.e. solutions which depend neither on the initial conditions nor 
on the precise form of the nonlinear term at large enough evolution time. In the example of the deterministic case 
(without noise), we display a sketch of the traveling wave solutions of the Id FKPP equation on Fig[TJ For sFKPP, 
the stochastic form of the 1-d equation, following a series of applications to QCI)_@, the recently found solutions 
(2^ can be interpreted as a stochastic superposition of traveling waves (see also (24 |) . 

Our aim is to look for similar properties for the 2-d version of the FKPP equation and find their consequences for the 
high-energy elastic amplitude, solution of the RFT. Our main results is the prediction of an asymptotic "universal" 
scaling form of the soft elastic amplitude T[Y, b) and in particular the prediction of an asymptotic expression for the 
expanding impact-paramter disk. 

The study of the general 2-d sFKPP equation ([8]) is interesting in itself and some results have been obtained already 
in the statistical physics literature (e.g. studies on the instabilities of the front wave [1^). In the following, due to the 
rotational symmetry in impact-parameter space of the elastic amplitude^ we shall concentrate our analysis of Eq.® 
to radial amplitudes, i.e. depending spatially only on the radial coordinate r = |f^. A comment is in order at this stage. 
Indeed, one could consider non azimuthally symmetric fluctuations contributing to a symmetric average. However, 
both the stochastic and the nonlinear character of the equation seems to invalidate this possibility. In particular 
anisotropic evolution may be caused by the noise [1^ and thus lead to a forbidden azimuthal symmetry breaking of 



^ Note that non azimuthally symmetric fluctuations could play a role in diffractive inelastic amplitudes, or other amplitudes which are 
not constrained by rotation symmetry. 
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the solution . We will assume that, if physical azimuthally assymetric fluctuations of the amplitude may exist, they 
are azimuthally averaged after a characteristic time much smaller than the typical evolution time in rapidity. 
One then obtain, after obvious integration over the azimuth, the following equation to be studied 

Uit,r)=drr U+-dr U+ U- {/^ + g , MlEEI „(t,r) , (9) 
at r V 2Trr 

where a curvature term ■^8^11 appears in addition to the original 1-dimensional FKPP equation Note also the 
modification of the noise strength by a factor l/27rr reflecting the symmetry constraint on the fluctuation strength. 
The white noise satisfies {iy{t, r), i'{t' ,r')) = 5{t' —t)5{r' —r). Looking for the asymptotic universal solutions of Eq.([5]) 
in terms of circular circular traveling waves in impact-parameter space is the goal of our paper. 



III. CIRCULAR TRAVELING WAVES: DETERMINISTIC CASE 



Let us first consider the equation ([9|) without the noise term, i.e. e = 0, namely the radial extension of the 2-d 
deterministic FKPP equation. It corresponds to neglecting the Pomeron loop contributions in the high-energy elastic 
amplitude. A series of results have been obtained for the 1-dimensional FKPP equation. The same methods, which 
we will adapt for the radial case, lead to new results. Indeed, the radial case is adding the derivative term ^drlJ to 
the standard FKPP equation and work in the half line r € [0, oo]. The deterministic equation to be studied is thus 

^ U {t,r) ^ drr U +- dr U + U - . (10) 
dt r 

We will prove the existence of traveling wave asymptotic solutions, which appear now as circular traveling waves 
(picturally reminiscent of those created by a stone falling in water). 

Let us recall first the main guiding principles of the FKPP traveling wave analysis. One distinguishes different 
regions, starting from the forward towards the backward of the wave, namely: the "very forward" region, the "leading 
edge" , the "wave interior" and the "saturation" regions. The universality properties mainly pertain to the "leading 
edge" region and its transition to the "wave interior" . In order to fulfill these universality conditions, characterizing 
the "critical" regime whence the traveling waves are formed, the initial conditions should be sharp enough in impact 
parameter, i.e. U{t = to,r^ tq) < e~^. This condition is fulfilled by considering an initial Gaussian form e^'' Z'*^ in 
impact-parameter^. We will keep this point-of-view in the following by considering e.g. a supercritical bare pomeron 
equipped with a Gaussian form in impact-parameter. 

Indeed, in the critical regime of high rapidity, the "very forward" region is driven by the initial condition, while 
the "leading edge" one develops a "universal" behavior [17| where three terms of the asymptotic expansion of the 
amplitude do not depend either on the details of the initial condition or on the nonlinear damping term. The "wave 
interior" posseses an exact scaling property (see further) while the "saturation" region depends on the nonlinear term. 
So, for completion, we will also give some hint on this region in the particular case of the initial RFT (k — 1) with a 
triple Pomeron coupling. 

In fact the main universal property of the traveling wave solutions in the deterministic case is the scaling property, 
namely 

Uit,r) = U{r-r,it)) , (11) 
where the "time" dependent radius rs{t) plays the role of the "saturation scale" in QCD Q. 



A. Universal "Leading edge" region 

Let us now derive the traveling wave properties in the radial case. For the leading edge, following a similar procedure 
for the 1-dimensional problem [3, [23| , one introduces in PU| an ansatz 

U{t,s^ r-vd) (X exp [-7,(5 + c{t))] T G ('-^^) , (12) 



It corresponds by Fourier transform to an exponential e T^i.e. a simple diffraction peak in transfer momentum for the elastic 
cross-section. 
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where Vc {resp. 7c) is the critical wave velocity {resp. critical slope) of the traveling wave front and c{t) describes the 
sub-asymptotic correction to the velocity v{t) = + dc{t)/dt. This ansatz describes the "velocity blocking" due to 
the critical mechanism. The point is that, being situated in the forward region where the non-linear terms in (fTO|) 
may be neglected, the form of the ansatz can be deduced from the linear part of the deterministic equation (jlOp . 
The only effect of the non-linearity is to ensure the "velocity blocking" by the compromise between the fast moving 
very- forward regime and the damping due to the non-linear unitarity bound (see, e.g. [§]). 

Inserting the ansatz in the equation (|10p and neglecting the small contribution from the nonlinear term to the 
leading edge, we can verify the equation for the dominant terms (successively in t*^, t~^) of the time expansion, 

see Appendix A. Note that the condition G(z) — > z when z ^ 0, is required in order to match with the scaling region 
(called the "wave interior" in [13 )• 

Adapting to the radial case the standard procedure [H, [13] , one finds 

U{r-rs,t) {r-rs) exp|-(r-rs)- '''' 



4 t 

r, = Vet + cit) ^ 2t - 2logt , (13) 

where is the average time position of the wave front (or "saturation scale" in the language of QCD M). Note 
that the form of the leading-edge front is the same as the one obtained in the 1-dimensional problem [ol.Tlq while 
the saturation scale is Vs = 2t — 21ogi instead of rg = 2t — |logi, due to the contribution of the new term ^drU 
characteristic of the 2-dimensionality of the initial physical picture. The saturation scale evolution is thus slower by a 
logarithrnic factor ^ \ogt. This 1/2 shift is due to the purely geometrical "curvature contribution" of the 2-dimensional 



problem [25| which combines with the coefficient 3/2 of the FKPP solutions. A third (and last) universal term in Vg 



behaving as t can also be derived and will add some new curvature contributions. 

The scaling (|lip is recovered from (IT^ in the region ^ ^ 1 giving rise to the simple expression 

U{r~rs) (r-rs) exp[-{r-rs)] , (14) 

where it ensures the transition with the "wave interior" domain. 



B. Circular wave properties for a triple Pomeron coupling 

1. Deep "Wave interior" region 

Taking into account now the specific quadratic nonlinear term of (jlOp (reflecting the original triple Pomeron 
coupling), we can explore the deep "wave interior" regime, adapting a method [l^] used for the similar problem in 
the QCD case Si]. 

Considering a scaling ansatz with an expansion in a small parameter to be determined by consistency with 
the scaling form p^ . 

U = U[z= ^— ^ = C/o + + C/4 + • • • , (15) 



A 



one obtains (see appendix B) 



with 



h A ^ ^ log ^ 

1 + (1 + e^)^ 4e^ 



U = + A-^ " . log Zl + ■ (16) 



..^^i^^^, (17) 

where the equality isobtained by matching^ at high enough t with the critical velocity of the leading-edge solution, 
namely A = 2. Note that it is easy to determine higher order terms by a system of nested linear differential equations. 



The matching between I I17I I and the leading edge velocity (1131 is not exact at sub-leading level, see Appendix bf B. 
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2. "Saturation" region 

The circular traveling waves being concentric around r = 0, the saturation region at small r is naturally expected to 
be different from the 1-dimensional one depicted in Fi^ where saturation starts from — oo. This is also made explicit 
by the drU/r term in Ea llOl which can no more be neglected or considered as giving second order effects as in the 
previous regions. For describing the saturation region, it is useful to go back to the full 2-dimensional form ([8]). It is 
convenient to introduce the S-Matrix element S — 1 ~ T, which is expected to be small in the saturation region. In 
those terms the deterministic^ 2-d equation writes 

S{t,r}=VlS- S+ , (18) 

Taking into account that is negligible, Eg. p^ boils down to a linear equation whose radial solution is easy to 
obtain if one notices that e^^S = is a solution of the two-dimensional heat equation, namely ^ W — W. A 
simple, azimuthal-invariant solution is thus: 

S{t,r) = l-U{t,r) = e-^^ W{t,r) = er^ {a-he-'^'^ , (19) 

where the constants a, b have to be determined by matching with the wave-interior region. This result shows the 
general feature of a time evolution towards the black disk limit 5 = 0, at large t. It does that in a non scaling way, 
since the approach to the black disk is not characterized by a single function r — ra{t). 



IV. CIRCULAR TRAVELING WAVE: STOCHASTIC CASE 
A. Quantum fluctuations and the Langevin equation 

As known from the seminal studies of Ref . [l6j , the effect of even very small fluctuations has an important impact on 
the solutions of the sFKPP equation. They may drastically modify the solutions of the sFKPP equations compared 
to the deterministic FKPP ones described in the previous section. Indeed, in the standard sFKPP case, one has 
been able to analyze [1^ that the small noise contribution has two superimposed effects. In the 1-d case, the typical 
expansion parameter appears to be not e itself but 1/loge, that is the inverse logarithm of the noise strength. At 
first order (starting in fact as 1/ log^ e) the correction has negative sign and corresponds to an effective cut-off on the 
amplitude as in [l^|. At the next order 1/ log'^ e, a positive contribution comes from rare but large fluctuations of the 
noise. 

In fact we shall now show that formula (jSOp for the analysis of the wave interior in the deterministic case the 
circular traveling waves with noise can be analyzed in a similar way than for the Id sFKPP case. However, some 
modifications will be due to the radial extension. Indeed, considering the initial Langevin equation Eq.Q, and the 
relation between the noise strength and the effective cut-off approximation [23l . [2^ , we are naturally led to an 
effective noise strength 

C(t) = e [2^r,(i)]-i/2 , (20) 

where in (fTU|) we have substituted —t in the expression of the cut-off. Indeed, this approximation can be 

justified by the accompanying factor U{1 — U) ^ outside r ^ rg. We see that for the radial case, the effective noise 
strength depends itself on the saturation scale and thus will possess a rapidity dependence. 

In fact the geometrical meaning of the noise strength ()20p is quite transparent. It takes into account the fluctuations 
at the periphery of the expanding disk in an azimuthally symmetric way. As an important consequence, the rapidity 
dependence of the effective noise will play an important physical role role, both at weak and strong noise regimes, as 
discussed now. 



In fact, the noise term would not play a big role anyway, since it is expected to have small effect in the "dense medium" characteristic 
of the saturated phase. 
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B. Stochastic traveling waves : weak noise 

Let us solve the weak noise regime of ([9|). Noting that choosing the variable z = r — J* dt'v{t') = r— 2t + ^ log(4t— 1) 
allows to take into account the radial term in the deterministic part of ©and to match the 2d radial case with the 
standard Id case (up to the modification ([20]) of the noise). Hence, taking into account the parallel properties of the 
radial equation with the Id, it is justified to export the detailed results obtained for the Id sFKPP equation [2^. 
However an important modification of the discussion for the radial configuration will appear due to the time-dependent 
noise strength (PO)) . 

The detailed effect of fluctuations has been derived 23] and leads to the following results: 
The result for the stochastic average over the amplitude is given [sil, [13] within some approximation (2§] by 



{U{r, t)) oc erfc 



DVt J ' " V 4 
where erfc(x) is the complementary error function and 



cxp — (r - rs) 



2 - erfc 



r-rs _ DVt \ 



D = 



31og=^(^) 



(22) 



(23) 



is the stochastic dispersion of the front. A complete description of the stochastic front, re-summing over all higher 
moments of the amplitude at weak noise, can be found in Ref.[29|. 

On a more general ground, as shown in [sil] and suggested by numerical simulations for the QCD case [11], one 
predicts a structure of "diffusive scaling" , namely 



U{r,t) 



(24) 



where the parameter I? is a characteristic diffusion coefficient, which may differ from the asymptotic (j23p . All in 
all, the solution of the stochastic equation can be understood as a dispersive distribution of event-by-event traveling 
waves with dispersion D. The random superposition of traveling waves transforms the "geometric scaling", valid for 
each of them into a "diffusive scaling" property (|24p for the average defining the final solution for the amplitude. 
However, following the numerical studies in the framework of QCD [32], diffusive scaling may require some evolution 
time to develop a sizable diffusion coefficient and thus may not be distinguished from "geometric scaling" at physical 
rapidities. 



C. Stochastic traveling waves : strong noise 

When the noise strength is tuned to increase, one observes a strong decrease of the average wave velocity, with a 
neat change of regime in the vicinity of a (normalized) noise strength of order one, see Fig. 2. Following Ref.[33|, the 
overall properties of the strong noise regime are as follows 

The solution of equation ^T0\ is a stochastic average of traveling waves at an average speed 



2 4TTr,it) 

c 



. ' (25) 



where ( is the normalized noise strength defined in ([^0]) . Hence the saturation scale follows from the equation 



(26) 
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noise strength 



FIG. 2: Average wave speed as a function of the noise for the sFKPP equation. Vertical axis: v/vc is the average traveling 
wave speed normalized to the speed Vc = 2 oi the deterministic FKPP equation; Horizontal axis: dimensionless noise strength; 
Dots: numerical results; Left Ime: weak noise analytic prediction; Right line: strong noise prediction. One observes (and may 
derive a maximal speed around a noise strength of order 10. The figure is from ref . [2C| . 



where the rapidity dependence of the radial noise plays the important role. Note that the limiting speed condition 
anyway requires v < Vc = 2 and thus from (j25p ATrrs{t) < e^. Hence the exponential behavior of (j26p must break 
down before the time evolution reaches the limit defined by 27r e^^*^'^ ^ e^. In fact, due to the rapidity decrease of 
the effective noise C of (|20p , the strong noise regime transforms progressively into the weak noise one, following from 
right to left the velocity curve depicted in Fig. 2. 

In the strong noise regime, there exists [33] an analytic solution for the average solution of the evolution equation 
(fTU|) . namely 

{U{t,r)) = i erfc^''"''^^*^ 



2 \ 2Vi 

1 (f-r (t)]^ 

' dre{r~r)exp-^- ±2^^ (27) 



2%/^7-oo 4t 

where erfc(a;) is the complementary error function. 

This result confirms the decrease of the velocity with increasing noise strength. It contrasts with the speed obtained 
in the weak noise limit by perturbative analysis around the F-KPP speed ~ 2 — 7r^| log(2C^)|~^. The expression ((27)) 
shows that the amplitude could be obtained from a superposition of step functions around = v t with a Gaussian 
form of width \j2Dt. The interesting point here lies in the dispersion coefficient: in the weak-noise analysis, it behaves 
like I log(2C^)|~'^. We have thus shown that the dispersion goes to a constant value 2 when the noise becomes strong. 
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V. THE POMERON AS A CIRCULAR TRAVELING WAVE 



Let us investigate the implications of the travehng wave properties on the soft Pomeron, within the framework 
that the soft interaction dynamics at high energies be governed by a "supercritical" bare Pomeron input. As we have 
shown in the previous theoretical sections, the circular traveling wave solutions are expected to appear due to the 
combined effect of the high energy evolution and of unitarity, which we will assume to be saturated^. Our analysis 
will be concerning the asymptotic regime of the circular traveling waves, leaving for further study the transition to 
this regime. 

We have seen that the evolution towards the saturation limit may depend on the non-linear terms, see Eq. p^ . 
Those terms may be physically more complicated than the single quadratic term of Eq. ([5]) . Hence we will focus on the 
"universal predictions" ,6.17. those which do not depend on the initial conditions and/or the structure of the non- linear 
damping. On the contrary, the parameter k, which is not fixed by the unitarity constraints, is the relevant parameter, 
playing an essential role in the Pomeron properties. 



A. "Phase diagram" as a function of noise 



The first step is to discuss which evolution regime we have as a function of k. For this sake, the noise strength ([^0]) 
can be conveniently written, restoring the Pomeron variables ([7]) 



2nr„ 



2n b. 



(28) 



It is important to note the following feature of the noise strength directly related to the (2-|-l)-d property of the RFT 
problem: it decreases together with the expansion of the impact-parameter disk and thus evolves towards weaker 
noise. However, this decrease, being governed by the evolution of the disk may be slow. 

The basic relation we will get comes from the structure of the wave speed reproduced in Fig. 2. It is obtained for the 
1-d case, but it happens to be indicative also for the radial case, whose universal properties are essentially similar, as 
we shall see. In Fig. 2, one may distinguish how the three different regimes we have analyzed in the previous sections, 
namely the zero, weak and strong noise respectively, can be identified on the plot where the K-dependent normalized 
speed v^/vc is displayed as a function of the normalized noise strength With our notations and using ([28)1 . we write 
by straightforward relations 



Vc 



1 



Vc 



~dt 



K 1 2dbs 
4^ C bsdY 



(29) 



where we have denoted v^ = ^ the actual wave front velocity and Vc — 2, the deterministic critical speed. Note 
that we have made use of (|28|) to substitute the bare Pomeron parameters \/ fia' by its expression in terms of the 
normalized noise. Our final expression thus writes 



Vk, 
Vc 



47r 



(30) 



where 6 = 2dbs/bsdY — d\ogA/dY, where A — nbl is the area of the effective impact-parameter disk for the collision. 
The obtained expression shows directly how the noise strength k parametrizes the normalized-speed vs. normalized- 
noise relation depicted in Fig. 2. It allows one to relate the "phase diagram" defining the different regimes of the radial 
sFKPP equation to a physical soft Pomeron feature, namely the exponent 5 of the expanding disk area A. 
When interpreting Fig. 2, one may distinguish the different regimes as follows using relation (|29)) : 



• The weak noise regime: 



C < 10" 



.9 < 



< 1. 



47r 



< 



IQ- 



dh^ 
dY 



(31) 



As we have seen in section III, the saturation limit at 6 = is not T = 1 at finite Y , contrary to the 1-d problem, see Fig[T] It reaches 
T = 1 when Y ^ 00. 
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• The strong noise regime: 

<>1.4 <-^5S^5 2V7^4»..C*>2§, (32, 

where we made use of the exact relation at strong noise To complete the picture, one adds 

• The zero noise regime: 



1 db. 



C «< 1 ~ 1 ^ ~ ^ -Sbs= ^ . (33) 



• The middle noise regime: 



.1<C<1.5 .2<^<.9 ^< — <^ < 2Va^ < 2— ^ . 34 

We see from relations pm34p that the parameter k plays the role of the order parameter of the RFT. Once given 
the physical observable (5, one knows the phase (c/. evolution regime) of the system from the determination of k. In 
particular for the original RFT action ^ , the phase is completely specified by S. Note also that the "bare" parameter 
2^0,' [I corresponding to the maximal critical speed of the disk relates to 6s<5/2 = 4y which is a "dressed" parameter 



dY 

in terms of a field theory. 

The question of determining the soft Pomeron properties thus boils down to the determination of b. We postpone 
a detaileed phenomenological study for the future"'^", but it is not too difficult to evaluate the order of magnitude 
of (5. Indeed, if the black disk limit would have been nearly reached, one would expect a geometrical cross-section 
a tot ^ A = irb^ and thus S ^ dlogutot/ dY ~ .08 where the last number is the well-known popular determination 
[lo| . This value for 5 should be considered as a maximum at present energies, since the black disk limit seems 
not to be fully reached (see, e.g., where one obtains smaller values of order (5 ~ (1 — 3)10~^). as an example 
in appendix C we show the phase diagram characteristics when choosing the conservative values of (5 ~ 10~^ and 
bs - 1 fermi = 5 GeV"\ 

In any case, one interesting remark is that for the whole range (5 S [1 — 8]10~^, the original RFT with k — 1 seats 
within the limit of the weak noise region. By comparison, the order parameter n takes a factor 100 in the interval 
between the weak and the strong noise regimes^"'^. Such high values of n are not it a priori forbidden, even if far 
from the original RFT action with a single triple Pomeron coupling. We will discuss in conclusion a possible QCD 
interpretation of these large values of n. 



B. Properties of the "wave front" 

1. An expanding impact-parameter disk 

Having now identified the phase diagram of the RFT solutions, we are able to discuss the characteristic features 
of the soft Pomeron as a circular traveling way by recasting the results of section II (resp. Ill) for the deterministic 
[resp. stochastic) traveling waves in terms of the physical variables through the relations ([7])- 

As a general result, valid in all cases, we find that the front of the traveling wave is situated around the impact- 
parameter value 

bs{Y) = ^^rs{t = ^lY) . (35) 

The expanding impact-parameter disk is thus related to the increasing function rs{t). This is the analogue of the 
rapidity-dependent "saturation scale" ^19.] discussed also in the framework of QCD traveling waves Q. Let us examine 



A first qualitative phenomenological exploration is discussed in Appendix D. 

It has been shown using field theory arguments l34ll that a maximal noise exists at f = Stt beyond which the traveling waves stop and 
the system does no more "percolate", i.e. the disk stops expanding with rapidity. 
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the equivalent "saturation scale" of the supercritical Pomeron. As we have seen previously, it depends on the phase 
diagram. Following Egs. p^ . (PT|) and ([^5]) respectively, we obtain 



bs{Y)-b,{Yo) = 2^/^iY-Yo) 
bs{Y) = 6,(^0) exp 



(y-yo)-^iog^ 



for zero noise 
37r2 log log (2C-2) 



An 



, for strong noise. 



, for weak noise 



(36) 



The middle-noise regime does not posses an analytic expression but its numerical implementation is possible and 
shown (in the reduced variables) on Fig. 2. The dots (• • • ) stand for sub-leading and/or non-universal higher order 
terms. Note that the constant terms implied by the initial condition at Y — Yq are also naturally not constrained 
by universality. The validity of relations p6p thus require a large enough interval Y — Yq. For all cases one finds 
from ([36]) that the disk expands with rapidity. For the "zero noise" and "weak noise" cases the asymptotic velocity 
is 2y/c/Jl, which is the critical velocity, and thus driven by the deterministic equation. The situation is different at 
"strong noise" where the solution is still moving linearly with rapidity Y but with a velocity governed by the order 
parameter k. 

For the "zero noise " and "weak noise" cases, we have thus identified new universal terms^^ not depending on the 
initial conditions nor on the specific non-linear damping terms. The existence of an universal rapidity expansion due 
to the supercritical Pomeron is, to our knowledge, a new result allowed by the Langevin formulation of the RFT and 
its traveling wave solutions. Moreover they appear to be quite different depending on non trivial phase diagram: in 
the deterministic case, the first (negative) correction to the radius oc Y behaves like log Y, while it is of order 
Y log^^ Y for "weak noise" , and thus a priori quite more important than in the deterministic case. 

For "strong noise", the obtained exponential behavior would not lead ultimately to a violation of the Froissart 
bound since the rise is tamed by the boundary of the strong noise regime. From (|28p and the relations ([32]), one has 



27r&, 



Vk/Vc < .5 



< 



47r 



(37) 



In fact this limit on the impact-parameter disk reflects the K-dependence of the noise strength and thus the evolution 
from strong noise towards weak noise through an intermediate middle-noise regime. Interestingly enough, this would 
mean for the radius (and thus for the cross-section near the black disk limit) a gradual transition from an exponential 
towards a squared logarithmic behavior in rapidity and thus an asymptotic restoration of the Froissart bound. 



2. Impact-parameter Scaling 

The scattering amplitude is related to the front profile of the dominant asymptotic traveling wave solution through 



T(Y, b)^U (t^ nY, r 



(38) 

at least in the region where universal results apply. U{t,r) being given by formulas PT|) . ([M)) and (P7)) respectively, 
one obtains 



T(Y, b) - T{b~ b^Y)} , for zero noise 
b-bs{Y) 



T{Y,b) - T- 



D^(Y)^a'{Y-Y,) 



, for weak noise 



T(Y, 6) - tI ^, -. \ , for strong noise , 

[Ds^a'iY^Y^) ' 



(39) 



(40) 



Note that a third universal term, behaving as 1/VY is expected to exist in the asymptotic expansion of the deterministic case from 1-d 
studies [gl llTll . We leave its determination in the present 2-d case for further study. 
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where the values of the front scales bs{Y) have to be chosen from (|36p for each corresponding individual regime. 
and Dg are the dispersion parameters for weak and strong noise respectively. We have introduced a rapidity scale 
Yi which denotes the effective rapidity where the dispersion of the noisy traveling waves becomes sizable. Indeed, 
numerical simulations for the QCD case 32] show that such a threshold do exist. Below that threshold the scaling is 
similar to the zero noise case. From (l23l) and (l27t. one finds 



27r2 2tt^ 

= = = 5- —— ; Ds = 2 . (41) 

31og3[2C-2] 3\og^ [AnbjK^^ ^ ' 

It is clear that the dispersion parameter D^^ cx 1/ log"^ Y may be quite small and thus one would then recover the 
same scaling as the zero-noise regime, but with the different weak- noise evolution bs(Y). Moreover the dispersion is 
proportional to ^a' which may be small even for a sizable value of \J a' \Ji for larger /i. 



3. Front profile 

In the regions where there exists a universal form of the wave front profile, i.e. within and forward to the wave 
front (the previously called "wave interior" and "leading edge" regions), one finds from formulas (fTTj) . ((22)) and ([27| 
respectively, 

T{Y, b) ^ {b- bs) exp | - " b^) - 4^^^^ (Y^_l^Ya) | ' ^^^^^ 



T(Y, b) - erfc I ^ i + exp 

\D^^a'{Y-Y,) j 



Jo f / b~b, D^^fi{Y-Y,) \\ . 
X < 2 — eric — > , tor weak noise 



D^y/a'{Y~Yi) 2 
2^a'{Y^) 



T(Y, b) ^ erfc < — ^ — > , for strong noise , (42) 



where D^, , Ds were given in ([4T|) . 

Some comments are in order about the front profiles. We may note that T{b, Y) for the zero noise case^^ goes to zero 
with b — bg, corresponding to the way how saturation is imposed as an "absorbing condition" [soj to the leading-edge 
approximations of the traveling waves' tails. Around and below b — bg the solutions get corrections to the spurious 
zero. However, these corrections (see, e.g. (|16p) are more dependent on the specific form of the nonlinear terms of the 
equation. 

Note also that the diffusive scaling (|24|) is exact for the strong noise case and with exact dispersion parameter 
Ds = 2. It comes from an exact solution of the statistical mechanic picture [1^. Indeed, the strong noise regime can 
be interpreted as the stochastic superposition of travehng waves of the simple form 9{b — bg), suggested by Eq. ip?)) . 
Moreover, at strong noise, the average solution is dominated by the strongest fluctuations, together with all correlators, 
as shown in [3311. 



VI. CONCLUSIONS 



We can summarize our results as follows; 

i) Reggeon Field Theory, when the bare Pomeron is "supercritical", can be formulated as a stochastic equation 
which is in the same "universality class" than the 2-dimensional version of the stochastic Fisher and Kolmogorov- 
Petrovsky-Piscounov equation. In this framework, "time" corresponds to a rapidity evolution Y — Yq and "space" to 
the 2-d impact-parameter 6 of the hadronic collision. 

a) Thanks to the mapping to the 2d-sFKPP equation, one is able to find the asymptotic and certain subasymptotic 
solutions of the RFT which were beyond reach of the purely field theoretical methods used in the past. These solutions 



It is also true for the noisy traveling wave solution I I21I I for weak noise. 
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possess an appealing "universality " property, which means that they do not depend either from the initial conditions 
or on the peculiar form of the non-linear terms ensuring the unitarity constraint on the elastic amplitude. 

Hi) To our knowledge, it is yet the only example of a supercritical Pomeron theory preserving an universal behavior 
of the amplitude. Usually, the factorisation property of a Pomeron as a Regge pole, on which relies the standard 
universality arguments (see, e.g. (lo|). is expected to be washed out by interactions for a RFT based on a supercritical 
Pomeron. This "universality" property is recovered in a very different way, since it comes from a dynamical mechanism 
based on the "critical" phenomenon associated to the formation of circular traveling waves. 

iv) The "universality class" property remains valid when the splitting and merging Pomeron vertices are of ratio 
K 7^ 1, which was the RFT value based on the unique triple Pomeron coupling. Indeed, k, which is a measure of the 
strength of the "Pomeron loop" contribution, plays the role of the "order parameter" in the phase diagram of the 
RFT. Hence the original RFT (with k = 1) lies in a specific phase of the diagram. A rough but realistic estimate 
(assuming RFT to be physically applicable) places the original RFT in the "weak noise" regime of sFKPP. However, 
thhis choice is not dictated by a theoretical constraint. The full phase diagram should allow a model-independent 
discussion of this degree of freedom, if compared directly with the phenomenology. 

v) More generally, depending on the dimensionless noise strength k the phase diagram is shown to lead to three 
specific phases corresponding respectively to zero noise {k = 0), weak noise {k — 0(1)) and strong noise (k — O(IOO)), 
(plus an intermediate middle noise one (0(1) ^ k < O(IOO)), for which explicit asymptotic regimes of solutions are 
obtained in the front and in the tail of the impact-parameter disk. Some other results, valid in the whole range are 
obtained in the case of triple-Pomeron coupling. 

There are intringuing theoretical^'' lessons to be drawn from our results The striking theoretical feature of the 
approach to the RFT through the mapping to the 2-dimensional sFKPP equation is its ability to avoid the complica- 
tions of an usual field theory formulation in the case of a "supercritical" bare Pomeron. Indeed, if at first a "critical" 
Pomeron theory for which the renormalization group exists raised some hope (see e.g. 0) it led to unphysical re- 
sults for hadronic reactions such as total cross-sections behaving as Y^^^"^. For a "supercritical Pomeron" the field 
theoretical methods, interesting as they may be (see e.g. Q), appeared to be technically complex with difficulties to 
conveniently handle the solutions. It thus seems that the traveling wave methods devclopped in the present study are 
well suitable for avoiding the obstacles. It is quite remarkable that, even a domain dominated by very large "quantum 
loop" contributions such as the strong noise phase discussed here, can be handled in a quite economic way. 

It is useful to list a series of interesting theoretical subjects which lie beyond the present study. One first problem is 
to get rid of the approximations made for the derivation of the solutions, the main one being to have replaced in the 
original equation ^ the r-dependent coefficients by their value on the wave front r^. It would improve the analysis 
to solve, even numerically, the original equation to check the validity of the approximation. 

The question of the azimuthal dependence of the noise is perhaps challenging and thus interesting. Indeed, It could a 
priori be possible that azimuthal symmetry be recovered after averaging over the noise. However, it is known from the 
study of a plane front ^2^] that instabilities may be created by the fluctuations beyond some threshold. It is reasonable 
to expect that the unitarity constraint should cut-off such inhomogeneities, leaving the purely radial solution valid. For 
non azimuthally symmetric variables, this is not so obvious and studies e.g. related to the "supercritical Pomeron" in 
diffraction dissociation or for particle production could lead to some interesting problems, such as the noisy structure 
of the diffraction disk [s^ . 

Finally, it would be natural and interesting to address the question of the relation of our results with QCD. A priori, 
there is a long way to go from an effective and thus "macroscopic" theory of Pomeron interactions to a "microscopic" 
point-of view based on quark and gluon interactions. Perhaps a tentative approach would be to notice that the RFT 
can be considered using the "hard" Pomeron as an input and thus depending of the value of a QCD coupling constant 
as. Since, finally, the only order parameter we have is n, this would mean that this parameter should be considered 
depending on as. It is interesting to note that the strong noise regime leads to the rapidity dependence bg oc e'*'^/". 
A "hard" pomeron behaviour which value of intercept is proportional to as would lead to choose k oc l/ag. Hence a 
"perturbative" property for QCD would be in relation with a highly quantum regime (large Pomeron loops strength 
k) in terms of RFT. By contrast a quasi-classical regime of RFT (zero or weak noise, small k) would be associated 
to a large effective coupling constant of order ~ 1/k. This speculative but intringuing "duality property" deserves 
certainly some interest. 



The present paper is not devoted to a phenomenologocal study. However, we have hsted in appendix D an outlook of possibly relevant 
phenomenological remarks. 
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Appendix A: Derivation of the radial "leading edge" 



In the following we shall restrict our analysis to radial amplitudes, i.e. depending only on the radial coordinate 
r = 1?^. One starts with 



^ U (t,r) = drr U + - dr U + U - C/2 . 
at r 

Let us introduce, following a similar procedure for the 1-dimensional problem [IB, the ansatz 

's + c{ty 



U{t,s = r-Vct) oc exp [-7c (s + cit))] T G 



(43) 



(44) 



where Vc {resp. 7c) are the critical wave velocity {resp. critical slope) of the traveling wave front. This ansatz is for 
describing the universal behavior of the wave in the leading edge region forward to the front [l7j . 

Inserting (|44|l in the equation ((43|) and neglecting the small contribution from the nonlinear term to the leading 
edge, we can verify the equation for the dominant terms of the time expansion, the different terms give: 



1 dU 
U H 
1 dsU 
U r 



7c - t 



G 
G 



G' 



7c " «c - Icc'it) + - + at-'^-^ [s + c{t)] — + 



G 



1 



S + Vet 



= \lc-t- 

(r_ G"^ 

~G^lP 



G 



Vnt 



(45) 



where the dots (...) indicate irrelevant sub-dominant contributions. 

Order by order in the late time expansion wc get the following relations 



= a7c/?c-a- — + ^ ^ (order t^^a) 
Vc G G 



IcVc = (7c) +1 
Vc = 2(7c) 
G' G" 



(order i°) 
(order t"") 



(46) 



with z = {s/t ") fixed and finite and c{t) = (3c logi. Hence we get the values of the critical parameters a = 1/2; 7c 
1; Wc = 2. The last equation of (|46p now reads 



= (r - 1) G{z) + - G'iz) + G"iz) 



(47) 



The condition G{z) z when z ^ 0, necessary to match with the scaling region (called the "wave interior" in [it'I), 
leads to /3c = 5/2 and G{z) oc z e~^ 1"^. One thus finally gets 



f/'-^-(r-r„t) ^ exp[-( 



r - r^)] (r - r^) exp I -- 



ti/2 



= 2t-21ogt 



(48) 



where I.e. stands for leading edge and Vs is the average moving position of the wave front (or "saturation scale" in the 
language of QCD [9]). Note that the form of the front is identical to the one obtained in the 1-dimensional problem 
d, [3 but the saturation scale is = s - c(i) + 1/2 = 2t - 21ogt (instead of = s - c(t) + 1/2 = 2t - ^ logi,) 
and thus slower by a logarithmic factor ^ logi. This shift can be interpreted (and checked) [25*1 as resulting from the 
superposition of the "curvature contribution" of the 2-dimensional problem with the dynamical slowing down of the 
Id FKPP solutions. 
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Appendix B: Derivation of the radial "wave interior" 

Let us introduce the formula (fTS]) into the equation pO|) and expand in powers of A^^, one gets 

1 



v{t) 

A Afdt'v{t') 
'v{t) 



U[ + Ull + (1 - 2C/o) C/i 



(49) 



where, one has replaced ^ drU 

In fact, we choose A = est. such that 



„t , ,, drU in ITUl) since the difference is Off) ^ At in the wave interior. 

J dt'vit') ^ * ^ ^ 



A/* dt'D(t') 



= 1, leading by simple integration to the equation 



rs = j = A<-ilog(Ar,)^ At-ilog(A2t-l) 



at large time. Solving the simple nonlinear equation of the first line of one easily gets 

1 1 



1 + e- 



1 + 



(50) 



(51) 



Knowing the solution for t/o, it is not too difficult to solve the linear equation, second line of ((49|) . obtaining with the 
appropriate boundary conditions 



U2 = o log ; 



(52) 



Appendix C: RFT phase diagram for ^ = 2 10"^ and bs ~ 5 GeV~^ 
• The zero noise regime: 



K < 1 



.025 GeV"^ 



(53) 



• The weak noise regime: 



C < 10"^ .9 < — < 1 k<2tt < .025 GeV" 



(54) 



• The strong noise regime: 



C > 1.4 (-^ = ^ < .5 K = 200 TT 

Vc 



The middle noise regime: 



.1 < C < 1.5 .2 < — < .9 

Vc 



/Jia' Sbs^ mbC > -05 GeV^^ , (55) 



27r < K < 200 TT .025 GeV"^ < ^/Jla' < .05 GeV"^ . (56) 



Appendix D: Phenomenological Remarks 

a) Scaling in impact-parameter. The main property of the traveling wave solution is its scaling structure in impact- 
parameter. On a phenomenological ground, assuming for simplicity a purely imaginary elastic amplitude, one obtains 
a scaling property of the elastic amplitude considered as a function of impact-parameter and energy. One may write 
ImTeiiY.b) ~ T(h — bg), where describes an expanding scattering disk A = nb'^ possessing universal slowering 
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corrections. Note that it is the "soft interaction" version, in the variables {Y,b), of a the similar "geometric scaling" 
[l^ of the "hard interaction" encountered in deep-inelastic scattering^^ and involving instead the {Y, logQ^) variables. 

As discussed in the paper, we expect scaling to stay approximately valid at weak noise, at least before a fully 
realized stochastic regime takes place where ImTeiiY. b) sa r|-^^^| . However, in this case, some non negligible 

corrections are expected to appear also in the disk radius, see p6p. 
h) Total cross sections 

The existence of an expanding disk in impact-parameter given by Eqs. (|36p may have a direct consequence on forward 
scattering amplitudes, and through unitarity, on the total cross-section at high energy. Indeed, taking as an example 
the ideal geometric relation atot A one finds for the zero noise case in appropriate units and for the dominant terms 
at asymptotic Y 

Y^-'-n2iL^...\ . (57) 
J 

Eq. (|57p saturates the Y'^ behavior given by the Froissart bound, and thus restores unitarity. However, our prediction 
is the existence of a "universal" YlogY correction term with strength governed by the bare Pomeron intercept /i. It 
is directly related to the correction term of the traveling wave speed 

Note that the result for the stochastic regime may be significantly different, namely 

a{Y) Y^ \l „^^- ^, > , for weak noise 

I W(^) J 

(j{Y) exp ^ , for strong noise , (58) 
n 

The results ((58|) call for comments. In the weak noise regime, it is clear that the stochastic corrections depending 
on the parameter k are of order 1 / log^ Y and thus significantly more important than the deterministic ones of order 
logF/F, see formula (1571) . Hence the Pomeron loop effects are expected, if their coupling k in the supercritical 
Pomeron scenario is effective, to have an observable effect. 

At strong noise, the behavior of the cross-section is entirely governed by the noise, with its characteristic parameter 
K. It is interesting to note that, if the phenomenological soft Pomeron of Ref.[l0| with intercept 1.08 is attributed 
to a strong noise scenario, it would correspond to k ~ IOOtt, which is discussed in the previous section UTTl Note 
that anyway the noise strength decreases like 1/^/Y and the strong noise regime will transform into the middle if not 
the weak noise one after some rapidity evolution. Hence the apparent violation of the Froissart bound will not be 
maintained at high enough energy. 

c) Modification of the large b behaviour 

It is interesting to study how the traveling wave behavior modifies the transfer momentum dependence of the elastic 
cross-section and thus the diffraction peak. Indeed, as we have seen previously, starting with an initial condition which 
is Gaussian in impact-parameter, the asymptotic traveling wave solution drives the solution of the evolution equation 
to a different, "universal" form {e.g. independent from the initial condition and the precise form of the non-linear 
damping terms in the equation). By Fourier transform, this evolution should change the structure of the amplitude 
in momentum transfer, and thus modify the diffraction peak. 

For an example we will start with the expression of the wave front in the "leading edge" domain at zero noise, 
formula . With some rearrangement of terms one can write 

6 + 2./^ iog(/i(r-ro))l' 

r(r,5)cx(6-6.) expAi(r-yo) exp-^ Aa'{Y~ Yo) ' ^^^^ 

In Formula (|59[) . we note that, apart the linear prefactor, the exponential behavior dominant at large Y boils down 
to the following modification 



Wriry|-exp|M>^-yo)-^^^^^ , (60) 



The term "geometric scaling" has been used long ago [35| for soft reactions at lower energies, but with a radically different formulation, 
namely T{b,Y) ~ T{b/R{Y)). 
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where 

6 = 6 + 2J-iog(/i(r-yo)). 

One recognize in the left-hand of (j60|) the solution of Eq.lH) reduced to the linear terms. Then ([60]) expresses the 
universal modification of the large b behavior of the amplitude due to the the universal leading-edge structure 

A similar study can be made for the stochastic case, e.g. from (|2ip. or even the explicit formulas (|22l27p . A detailed 
phenomenological study of all these aspects is deserved, based on the results of the present theoretical paper. 
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